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' Using a chirally symmetric Lagrangian, which contains quarks as ele- 

Qh! mentary degrees of freedom and mesons as bound states, we investigate 

the expansion and hadronization of a fireball, which initially contains 
' only quarks and produces mesons by collisions. For this model, we 

study the time scales of expansion and thermal and chemical equili- 
bration. We find that the expansion progresses relatively fast, leaving 
not necessarily enough time to establish thermal and chemical equilib- 
rium. Mesons are produced in the bulk of the fireball rather than at a 
surface, at a temperature below the Mott temperature. Initial density 
fluctuations become amplified during the expansion. These observa- 
tions challenge the applicability of hydro dynamical approaches to the 
expansion of a quark-gluon plasma. 
PACS numbers: 25.75.-q, 12.28.Mh, 12.39.Fe, 24.10.Lx 
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I. INTRODUCTION 



The physics of the quark-gluon plasma (QGP) has formed one of the most impor- 
tant points of interest in nuclear and high energy physics during the last years 
While possible evidence for its creation has already been obtained in Pb-Pb colli- 
sions at the CERN SPS, new experiments at RHIC and LHC are under construction. 
It is expected that these experiments will reproduce the SPS results and beyond that 
confirm other signals for the creation of a QGP, which have been proposed in the 
literature. 

With these improved experiments coming up soon, there is a clear need for a 
theoretical description, which is able to describe the different stages of a heavy ion 
collision, as there are (i) the formation of the fireball, (ii) its expansion within the 
QGP phase, (iii) the phase transition to the hadronic phase, (iv) the expansion 
within the hadronic phase, and (v) the decay of the fireball into noninteracting 
fragments. With the exception of the first of these stages, this is frequently done 
using hydrodynamical models due to the lack of a more realistic approach. This 
approach assumes, that one has local thermal equilibrium during the whole evolution 
at least for the light and thus dominant components of the system, so that it is 
sufficient to compute the space-time dependence of the energy-momentum tensor 
using a phenomenological equation of state, which may be extracted e. g. from 
lattice simulations. One of the drawbacks of this method is that it is implicitly 
assumed that the time scale for collisions is short compared to the expansion time 
scale, which is by no means obvious due to the large expansion velocity of the fireball. 
Furthermore, the equation of state is unknown for /i = 0, since this sector is still 
not accessible for lattice computations. Another weak point is that hydrodynamics 
has no intrinsic freezeout mechanism, so that the breakup stage into hadrons of the 
expansion must be described by additional assumptions. 

Another common approach for the modelization of heavy ion collisions are trans- 
port calculations based on cascade codes. These models have the advantage of being 
generic non-equilibrium calculations. They suffer, however, from the fact that usu- 
ally the in-medium modifications of particle properties and interaction cross sections 
are neglected. Furthermore, no attempt is made to apply these calculations to phase 
transitions. This limits the validity of these models to the late stages of the expan- 
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sion, where the system is already in the hadronic phase. 

The approach we present here differs from those presented above in that we 
attempt to construct a cascade type model, which is nevertheless able to describe a 
transition from a pure quark phase to a hadronic phase, and which also takes into 
account the in-medium modifications of quark and hadron properties. Doing this 
in a fashion which describes all of the phenomena is nevertheless impossible, since 
one of the major effects of the QGP transition, the confinement, is presently not 
understood sufficiently to implement it in a non-equilibrium model. It is, however, 
known that confinement does not play an important role in the low energy hadron 
phenomenology. This sector is rather dominated by chiral symmetry, which can 
be well described by effective models. Since the importance of chiral symmetry 
is well established and the behaviour of effective chiral models is also extensively 
studied, an incorporation of this symmetry into the dynamical theory of heavy ion 
collisions seems mandatory. We thus take one of these effective models, the Nambu- 
Jona-Lasinio (NJL) model Hf|, for which a non-equilibrium treatment has been 
recently derived and for which first numerical solutions have been reported 



in Refs. |p], [T0| , [TT]] . The advantages of this approach are obvious: Since the NJL 
model starts from a Lagrangian, which contains only quarks as elementary degrees 
of freedom and treats pions as bound states, one is able to describe the transition 
of a hot system containing only quarks initially to a system which contains quarks 
and hadrons. The technical tools for the treatment of this transition have been 
developed in Ref. ||. The possibility of including both quarks and hadrons, the 
former appearing in the initial state while the latter will be present in the final 
state, distinguishes the NJL model from most of the other chiral models such as 



chiral perturbation theory or the a-model ||12|| , while the description of hadronic 
properties at low temperatures is equally good for all of them. A comparison of the 
NJL mesonic mass spectrum with the lattice calculations of Ref. |13| shows that 
also at high temperatures the NJL model provides a good description of mesonic 
properties. Besides that, another advantage of our approach is that it is possible 
to derive the whole expansion scenario from one single Lagrange density. Since the 
NJL model does, however, not include confinement, a price has to be paid for this, in 
that the transition between the two phases is not given by a confinement transition, 



but rather by a Mott effect |L4]. As a consequence, one has thus free quarks at all 
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temperatures. 

The numerical method we choose for the solution of the transport equations is an 
algorithm of the QMD type ||15|| . In this approach, it is assumed that the particles 
can be described by (nonrelativistic) gaussian wave functions, which contain two 
time dependent parameters, the position of the centroids in coordinate and momen- 
tum space. The total wavefunction of the system is assumed to be the product of the 
wave function of the individual particles. The form of the wavefunction remains time 
independent and the time evolution of the parameters is obtained applying a time 
dependent variational principle, which yields time evolution equations very close to 
the classical time evolution equations. Being a n-body theory, this approach allows 
to investigate in detail the time evolutions of correlations and fluctuations, which 
is beyond the range of applicability of BUU type transport theories. For details we 
refer the reader to Ref. [|T6| . 

Previous approaches to a nonequilibrium treatment of the NJL model have been 
given in Refs. [H ITOlPl - This work goes beyond Refs. HT0| , [rT| in that these papers 
only considered a system of quarks moving within a mean field, whereas we consider 
also collisions. Reference || also included collisions in a relaxation time approach, 
however, mesonic degrees of freedom and the mechanisms of their creation were not 
studied there. This is thus the first work, in which the chiral phase transition is 
studied including collisions between quarks as well as a hadronization of quarks into 
mesons. 

This paper is organized as follows: In Sec. O, we describe our interaction model 
and its numerical treatment. Numerical results are given in Sec. [TT1]. We summarize 



and conclude in Sec. [V 



II. DESCRIPTION OF THE MODEL AND THE ALGORITHM 



The model which we use throughout this paper is the Nambu-Jona-Lasinio (NJL) 
model || in its two flavor version, which is defined by the Lagrangian 

,2 / - \2" 



C = ip (i — m ) vp + G 



(2.1) 



Here, ip denotes the quark fields, which are implicitly understood to carry flavor and 
color indices, and f are the Pauli matrices in flavor space. A current quark mass 
m is introduced, which provides a small explicit chiral symmetry breaking. 
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For a review of the equilibrium properties of this model, the reader is referred 
to Ref. ||. Here we summarize only those points, which are essential for our calcu- 
lations. At small temperatures and densities, the interaction in Eq. fl2.1|) leads to a 
spontaneous breakdown of chiral symmetry, which is described by a finite effective 
quark mass m q . To lowest order in an expansion in the inverse number of colors, 
1/N C ||17|| , which for quarks is commensurate with the Hartree approximation, m q is 
given by the gap equation 

m q = m + 2Gm q [ f P ) (2N c N f - n q (p) - n q (p)) , (2.2) 
J\p\<a {ivy tiq{P) 

where E q {p) = ^jp 2 + m 2 is the quasiparticle energy of the quarks and n q (p), n q (p) 
denote the quark and antiquark momentum distributions, respectively. Since the 
NJL model is non-renormalizable, we have restricted the integral in Eq. ( |2.2| ) to 
momenta smaller than a cutoff A. 

Pions appear in the NJL model as bound states of quarks and antiquarks in the 
pseudoscalar interaction channel. As a consequence of the Goldstone theorem, they 
become massless in the chiral limit, itlq — > 0. If, however, mo is finite, they also 
acquire a finite mass. Their quasiparticle energy can be extracted from the poles of 
their effective retarded propagator, A^(p ,p), which is given by [Q 

2G 

Agfa,.*)- T^onM • (2 ' 3) 

where Hp S (po,p) is the irreducible retarded pseudoscalar polarization function. The 
quasiparticle energy of the pion is thus given as the solution of 



l-2GU^ s (E n (p),p) = . (2.4) 



In the following we will make the approximation E n (p) = \Jp 1 + m^, where is 
computed as E n (0) from Eq. ( |2.4j ). The retarded polarization is computed, again to 
lowest order in 1/N C , from the diagram given in Fig. |I| In this case, the irreducible 
pseudoscalar polarization is explicitly given by 

n£ 5 (p ,p)= (2.5) 
d 3 q 2N c N f -n q {q)-n q {q) pq(p 2 + 2pq) - 2p 2 Q E q (q) 2 



1<a(2tt)3 Eq ffl (p 2 + 2pq) 2 -4p 2 E q (q) 2 + iep 

with the term iepo in the denominator taking care of the causal structure. 
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Due to the scalar coupling channel in the Lagrangian ( |2.1|) , the model contains 
also a scalar resonance, the a meson. Its mass can be computed from a dispersion 
relation similar to Eq. (|2.4|), replacing the pseudoscalar polarization by the scalar 



nf( Po ,p)= (2-6) 

d 3 q 2N c N f - n q (q) - n q (q) (pq + 2m\){p 2 + 2pq) - 2p 2 E q {q) 2 



/ b1 <A(27r)3 Eq ffl (p2 + 2pq) 2 -4p 2 E q (q) 2 + iep 

This function is computed from the same graph as given in Fig. [TJ, with ry 5 r replaced 
by one. 

The mass spectrum computed from Eqs. ( |2.2|) , (|2.4j) in thermal equilibrium is 
shown in Fig. |2[ The parameters used for this plot are mo = 4MeV, GA 2 = 1.989 
and A = 820 MeV. At zero temperature, chiral symmetry is spontaneously broken 
and the quarks appear as constituent quarks with masses of m q = 323 MeV. As a 
consequence, light Goldstone pions appear as bound states, which are, in the chiral 
limit, massless, but for the parameters given have the mass = 136 MeV. The a 
meson has a mass of m a = 653 MeV. As the temperature rises, the Pauli blocking 
term of Eq. ( |2.2j ) becomes more and more important. Thus the quark mass drops 
at high temperatures and chiral symmetry becomes restored. Concomitantly, the 
pion mass will rise, so that at the Mott temperature Tm its mass becomes equal to 
that of its constituents, tr^Tm) = 2m q {TM)- For our parameter set, this transition 
occurs at T M = 219 MeV. At temperatures higher than T M , the ground state of 
the model is no longer given by the pions, but rather by the quarks, while the 
pions themselves become unstable resonances. This instability models, to a certain 
extent, the deconfinement transition of QCD. For a more detailed discussion of 



this Mott transition the reader is referred to Ref. [TJj]. As a consequence of chiral 
symmetry restauration, the a mass becomes degenerate with the pion mass at high 
temperatures. 

It has been detailled in Refs. |^-|7|U, how the treatment of the NJL model is 
generalized to non-equilibrium. In particular, it has been shown, that the quasipar- 
ticle energy of the quarks is given by an equation similar to Eq. ( j2.2| ), where m q and 
the momentum distributions become space-time dependent, whereas the meson self 
energy in turn is given by to Eq. (|2.4| ) with a space-time dependent polarization. 
The equations of motion for the one particle distribution functions is then given by 
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[d t + d p E q d x - d x E q dp) n q (t, x , p) = I co ii, q [n q , n n } (2.7a) 
d t + d p E n d x - d x E n d p ) n n (t, x,p) 

-^coll,7r \P<qi Tin ] , (2.7b) 



where I co ii q [n q , n n ] and I C oii,n[n q ,n n } are collision integrals of the Boltzmann type. 

In the present paper, we solve Eqs. ( |2.7| ) using an QMD-like approach. To this 
end, we replace the one particle distribution function n q by the parametrization 

n q (t,x,p) = £exp (- [f ~ff )]2 ) exp (-^ Ip-pM?) (2-8) 

and analogously the antiquark distribution function n q and the pion distribution 
function n n . The normalization in Eq. ( |2.8|) is chosen in such a way that 

d 3 xd 3 p , t . . . 

-^n g (t,x,p)=N q , (2.9) 

where N q is the total number of quarks. The center points Xi(t), pi(t) in Eq. ( |2.8| ) 
are moving on the characteristics of Eq. ( [2.7| ) : 

—x i {t) = B p E (2.10a) 
dt 

— Pi(t) = —d x E + collision contributions . (2.10b) 
dt 



The second term on the right hand side of Eq. ( [2 . 1 (J b ) serves to describe the effects 



of the collision integral. It is computed as follows: for each pair of particles, the 
coordinates and momenta are boosted into the respective two particle center of mass 
frame. Here, the distance vector of the particles is decomposed into its longitudinal 
part, i.e. the part parallel to the relative momentum, and its transverse part, i.e. 
the part perpendicular to the relative momentum. We decide, which scattering 
processes are possible for the pair in question and compute the respective cross 
sections as well as the total cross section cr to t = Y, cr k- A collision happens, if (i) 
the longitudinal distance of the incoming particles in the CM frame is smaller than 
the distance traveled by the particle within the time step and (ii) the transverse 



distance is smaller than Ycr tot /7r. In this case, the specific process is chosen with 
probability cr k /cr tot . Neglecting the dependence of the differential cross section on 
the angles, we choose the direction of the outgoing momenta randomly in the center 
of mass system and boost them back into the original reference frame. 
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The processes implemented so far which are treated in this way belong to two 
classes. The first of these are elastic scattering processes of the form qq qq, 
qq <->• qq and qq <-> qq 0]. The generic Feynman diagrams of these processes in 
leading order of the 1/N C expansion are shown in Fig. |3]. One has two interaction 
channels, which both proceed via the exchange of a scalar or pseudoscalar meson. 
The typical cross section for these processes are of the order of millibarns. The 



second class is represented by hadronization processes, qq <-> nn [fL9| , P0|| , f° r which 
the Feynman diagrams are shown in Fig. [|. Here, one has an s-channel, which 
proceeds via the exchange of a scalar resonance, and a t- and u-channel. Also here, 
the typical cross section is of the order of millibarns. For more details and plots of 
the cross sections, the reader is referred to Refs. |JT8| |20|] . 

A third class of processes, which cannot be treated by the scheme detailled above, 
is the decay of pions into two quarks, for which the Feynman diagram is shown in 
Fig. [5]. This process is possible if the pion mass fulfills the Mott condition > 2m q 
||. Here we proceed as follows: first we compute the mean lifetime r of the pion 
due to the process 7r —> qq. Then we decide with probability 1 — exp(— At/r), where 
At is the time step of the calculation, if the particle decays. In this case we again 
choose the direction of the outgoing particles randomly in the meson rest frame and 
boost to the original frame. 

In general, the quark and meson quasiparticle energies as well as the scattering 
cross sections appearing in the collision integrals are complicated functionals of the 
particle distributions, as can be seen from Eqs. ( |2.2| ), ( |2.4j ) and the expressions for 
the cross sections in Refs. []T^-CT. An exact computation of these quantities in a 
non-equilibrium situation is thus a tremendous numerical task. Instead of making 
an exact computation at each time step, we thus take a shortcut by defining an 
effective temperature. This is done as follows: First, the gap equation ( |2.2| ) with 
the parametrization ( |2.8|) inserted is solved exactly at each particle position to give 
the quark mass field m q (x, t). Then the effective temperature is computed from the 
condition that the quark mass at the respective point is equal to the equilibrium 
quark mass computed from the effective temperature: 

m q (x,t) =m^(T eff (x,t)) . (2.11) 

From this effective temperature we compute afterwards the pion mass and the scat- 
tering cross sections, using the equilibrium expressions. Note that T eS is not the 
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thermodynamic temperature, but rather an auxiliary quantity for the computation 
of meson masses and cross sections. It becomes, however, equal to the thermody- 
namic temperature if thermal equilibrium is reached. 



Note that the second factor on the right hand side of Eq. (|2.12j ) is negative, so that 
the effective force acting on pions due to the mean field is opposite to the one acting 
on quarks. With these prescriptions at hand, Eqs. ( |2.10|) are solved to give the time 
evolution of Xi(t), Pi(t) and thus the time evolution of the one particle distributions. 

The initial condition is chosen in the following way: In the beginning, no pions are 
present, whereas the number of quarks and the number of antiquarks are equal. The 
spatial center points of the Gaussians representing the quarks Xi(0) are distributed 
homogeneously within a sphere of radius r^. The momenta Pi{0), on the other hand, 
are distributed according to a Fermi distribution at a given temperature T , which is 
a free parameter, and chemical potential /i = 0. The total number of quarks in the 
initial state is given by the momentum integral over the Fermi distribution times the 
volume. Since the NJL model is a low energy theory, the momentum distribution is 
cut off at p = A. 



In this section we show the numerical results of our simulation program. All 
computations have been performed using the parameter set mo = 4MeV, GA 2 = 
1.989 and A = 820 MeV. The width of the Gaussians in Eq. ( |2.8| ) was chosen to be 
w = 2fm. 



To give an overview over the time evolution of the fireball, we first describe one 
specific example of a system with initial radius r = 7fm and initial temperature 



A similar procedure is used for the mass gradients, which appear in Eq. ( |2.7| ). 
First, we compute dm q from an exact differentiation of Eqs. ( |2.2j ), fl2.8p . The pion 
mass gradient is then computed from 




(2.12) 



III. NUMERICAL RESULTS 



A. Evolution of the Fireball 
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T Q = 280 MeV. A snapshot of the expansion is shown in Fig. || The upper part 
of Fig. ^| shows the initial state, which consists of a sphere filled up with quarks, 
which are denoted by dark balls. At this time, the system contains 4684 quarks 
and antiquarks and has a kinetic energy density of about 1.8 GeV/fm 3 . The bottom 
part of Fig. |^ shows the same system 25fm/c later. At this time, approximately 
60% of the quarks have been converted into pions, which are denoted by the light 
balls. Since the NJL model is not confining, the remaining 40% light quarks will 
not hadronize and remain present in the final state. 

Figure [7| gives the time behaviour of the constituent quark masses, averaged 
over all solid angles, as a function of r at times t — 0, 5, 10, 15, 20 and 25fm/c. 
At t = 0, all quarks are sitting within a sphere of radius 7fm. In the center of this 
sphere, the particle density is maximal and the quark mass is thus low. Towards 
the surface, the density drops due to the gaussian shape of the distribution function 
and concomitantly the constituent quark mass rises. Note that the quark mass is 
only known at the centre points Xi(t) of the quark distribution functions and the 
plot ends thus at r = 7fm. For larger values of r, it would continuously rise towards 
the vacuum value, thus forming a sort of 'potential well', as was previously shown in 
Ref. ||. At finite times, the system starts to expand due to the thermal motion of 
the quarks. This expansion first depletes the surface of the fireball, while the particle 
density in the centre stays high, as can be seen from the plot at t = 5fm/c of Fig. |7|. 
At this time, some quarks have already gained their vacuum mass, while the mass of 
those quarks, which are in the centre, is approximately unchanged. At later times, 
when the depletion reaches the centre, the 'potential well' begins to flatten, until 
finally all quarks have the vacuum mass. At this time, the interaction due to the 
mean field ceases. This behaviour has been observed previously in Refs. @,|I0]. 

The effective temperature, averaged over all solid angles, is shown in Fig. 
Since this quantity is coupled to the constituent quark mass via Eq. ( |2.11| ), it shows 
qualitatively the same behaviour as Fig. [?]. Initially, one has a high temperature 
region in the centre of the fireball and zero temperature outside. As the expansion 
progresses, the temperature drops in the centre. At t = 7fm/c, the temperature is 
lower than the Mott temperature at all space points. At t = 25fm/c, one arrives at 
a temperature of approximately 50 MeV, which is sufficiently close to zero to give 
no mean field contribution any longer. 
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The result which can be drawn from this behaviour is that the high energy density 
in the initial state leads to a large expansion velocity, which in turn results in a short 
lifetime of the plasma. Although the increase of the quark mass should make the 
expansion slower, this does not make a large effect. Quantitatively, this can be 
seen from Fig. |9], where the mean radius is shown for three systems with identical 
initial conditions, but different expansion mechanisms: one with full interaction, 
one expanding according to a Vlasov equation (i. e. without collisions) and one with 
no interaction at all. The curves for the interacting systems lie below the one for 
the interaction free system, thus showing that the growing quark mass leads to a 
slowing down of the expansion, which is, however, not very strong. In the case of 
a fully interacting system, some of the quarks are converted into pions, which are 
not slowed down, so that in this case the mean radius is somewhat larger than in 
the Vlasov case. One may ask the question whether this short lifetime is a generic 
result or only a consequence of the approximation of the NJL approach. In view of 
the high initial energy density and the fact that the quark masses are rather small, 
it is hard to see how an expansion well below the speed of light can be achieved. 

Figure [TO] shows the time behaviour of the angular averaged particle densities, 
both for quarks and for pions. In the initial state at t — 0, only quarks are present. 
At later times, pions are produced and thus contribute also to the particle density. 

The time dependence of the particle multiplicities is shown in Fig. |Tl|. Here it can 
be seen that the production of mesons starts immediately after the beginning of the 
evolution. The pion multiplicity rises steeply at t = 0, while the quark multiplicity 
goes down. This steep rise can be understood from the observation, that the system 
is rather dense initially, thus giving rise to a high collision rate. At later times, the 
multiplicity curves flatten. There are two possible reasons for this flattening: either 
the back reaction txtx — > qq might become important due to the large abundance 
of pions or the particle density becomes to low to cause a further change of the 



multiplicities via collisions. To answer this question, we show in Fig. [12| the number 
of collisions per unit time, dN c /dt, as a function of time. It can be seen from this 
figure, that the vast majority of collisions happen before t = lOfm/c. If the reaction 
Tin — > qq became important, the number of collisions would stay high, while the 
multiplicities do not change any longer. Thus the reason for the flattening of the 
multiplicity curves is the breakup of the fireball rather than the approach to chemical 
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equilibrium. 



B. Meson Production Mechanisms 

In the last subsection it has been pointed out, that approximately half of the 
quarks in the system do not hadronize and remain present in the final state. In order 
to investigate this behaviour further, we have simulated systems with initial radii of 
3, 4, 5, 6 and 7fm. For each initial radius, two runs have been performed. The final 
multiplicities for each system, averaged over the two runs, are given in Tab. |. It can 
be seen from this table, that the final state for all these runs contains about 50% of 
quarks, but it can also be seen, that there is a systematic trend towards a higher pion 
abundance for larger systems. Whereas one has 42% pions at r = 3fm, the pion 
abundance rises to 60% at r = 7fm. This can be interpreted as an approach towards 
chemical equilibrium, which is nevertheless not reached, since in an equilibrated cold 
system, only pions should be present. The production of pions from quarks, however, 
takes a finite time. The lifetime of a small expanding fireball is too short to hadronize 
all quarks into pions without involving confinement. Nevertheless, the lifetime rises 
with the initial radius, so that more quarks can hadronize for a large system. 

The hydrodynamical picture of the plasma expansion relies on the assumption 
of local thermal equilibrium, which is the limit of infinitely large cross sections of 
transport theory. In this case, the transition from quark matter to hadronic matter 
takes place immediately after the temperature drops below the critical temperature. 
In an expanding finite system, however, one has to take into account that the finite 
size of the cross section leads to a finite mean free path and the hadronization takes 
a finite time. To verify whether in our nonequilibrium, finite particle number cal- 
culation a comparably strong relation between meson production and temperature 
as in hydrodynamics is realized, we have plotted in Fig. ^ the creation probability 
of pions versus the effective temperature. It can de seen that pions are predomi- 
nantly created at temperatures well below the Mott temperature within a range of 



approximately 50 MeV. This figure should be compared to Fig. 23 of Ref. [20|, where 



it has been shown that the mean time for the conversion of quarks into mesons is 



minimal within the same temperature range where Fig. 13 shows a maximum of the 



pion production. This means that the mean free path is minimal below the Mott 
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temperature and quarks tend to hadronize here. We thus conclude that the finite 
creation time of mesons cannot be neglected against the expansion time scales of 
the plasma, i.e. effects of the finite mean free path become important. 

This behaviour is consistent with the density curves of Fig. [R]. Generally one 
observes, that in the early stages, when the system is still hot in the center, the meson 
density is maximal at a finite value of r. At later times, when the temperature in 
the centre drops below the Mott temperature, one has the situation, that the meson 
density follows the quark density. 



C. Density Fluctuations 



The initial conditions used above are, besides statistical fluctuations in the distri- 



bution of the particles, spherically symmetric. We see, however, in Fig. |T|, where we 
plot the summed particle densities for quarks and pions along the coordinate axes, 
both for the initial and final state, that these small fluctuations become amplified 
during the expansion and that the system does not at all keep the symmetry. This 
is understandable because the increase of the parton mass with decreasing density 
is equivalent to an attractive force between them. Therefore partons tend to cluster 
and the initial fluctuations are the seeds for this clusterization. This observation 
challenges, like the finite mean free path effects, the possibility to employ hydrody- 
namical approaches, which are only applicable if initial fluctuations are damped. 



D. Variation of Initial Conditions 

1. Cylindrical Initial Conditions 

The initial conditions given in Sec. |I| are not the ones one would expect in a 
heavy ion collision. To come closer to the experimental situation, we made one 
run with cylindrical initial conditions. The spatial center points in this run are 
equally distributed within a cylinder which, in longitudinal direction, extends from 
-z max < z < z max with 2 max = 2fm and in transverse direction has a radius of 6fm. 
The distribution of the momenta perpendicular to the cylinder axis was chosen to 
be gaussian 
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p(p^) ~ exp (^-^ (3.1) 

with width 5p = 250 MeV. The distribution along the cylinder axis, however, was 
chosen to depend linearly on z in such a fashion, that p z = ±A is reached at the 
cylinder faces at z = ±z max : 

p(p z ) = 5{p z - Az/z max ) . (3.2) 

The total number of particles was 2000, which gives approximately the same particle 
density in the initial state as for the other runs. After 30fm/c, 1126 quarks have 
been converted into mesons, which corresponds to 56% of the total particle number. 



In Fig. [15], we show the time dependence of the multiplicity. This figure shows the 



same trend as Fig. |TT]. The geometry does thus not have a large impact on the 
particle multiplicities. Figure [16| shows the distribution of the center points Pi(t) in 
rapidity both for the initial and the final state. At t — 0, one sees a distribution 
with two maxima at y = ±1.1. At t = 30fm/c, the quark rapidity distribution has 
become flatter and narrower, while the meson distribution is flat at midrapidity, but 
has a larger width than the quark distribution. This comes from the fact that quarks 
are slowed down during to the expansion and thus move to lower rapidities. Mesons, 
on the other hand, have a lower mass than quarks at the time of their creation and 
thus obtain a higher momentum as the incoming quarks. This leads to a spread of 
their rapidity distribution. 

Figure [17| shows the mean transverse momentum as a function of rapidity. Here 
one encounters the same behaviour as for the particle density, in that the distribution 
for the mesons is wider than for the quarks. The geometry of the expansion does 
thus not influence the px spectra, but shifts the quarks to midrapidity. It is not very 
astonishing that mesons have the same mean transverse momentum in midrapidity, 
because the process qq — > nn leads to about the same value of px- 



2. Mixed Temperature Initial Condition 

In view of the short lifetime of the fireball, the question arises whether the lifetime 
is long enough in order to establish thermal equilibrium. For this reason we have 
studied a spherically symmetric system, in which half of the particles with r > Tq/2 
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obtained a momentum according to a temperature of 180 MeV ('cold' particles), 
whereas the other particle momenta were initialized with a temperature of 280 MeV 
('hot' particles). If thermal equilibrium was established, the 'cold' particles would 
obtain the same energy as the 'hot' particles in the course of the expansion. The 
result is shown in Fig. [TBI, where we plot the energy per particle for both subsystems 
as a function of time. It can be seen, that the energy of the 'cold' particles approaches 
that of the 'hot' particles, but stays below. One has thus no complete thermal 
equilibration. This reflects the already discussed observation that the expansion 
velocity is very large, too large to produce a global equilibrium. 



E. Variation of Hadronization Cross Sections 



Since the NJL model is nonconfining, one is faced with the question, if a suffi- 
ciently enlarged hadronization cross section could lead to a system, which complete 
hadronizes into mesons. We thus took the initial conditions of Sec. |Tl ID 1| and made 
several calculations, in which we increased all cross sections of the hadronization 
type, i. e. the ones for qq — > tttt as well as those for the back reaction tttt —> qq, 
artificially by a constant factor. The result is given in Table |TJ, where we give the 
particle abundancies in the final state for these runs. One sees that the increased 
hadronization cross section leads to an increased pion abundance in the final state, 
but nevertheless not to a complete hadronization. 



IV. SUMMARY AND CONCLUSIONS 



In this paper, we have studied the expansion and hadronization of systems, 
which interact according to a chirally symmetric Lagrangian. These calculations 
cannot serve to reproduce experimental data, since the underlying interaction is not 
confining. Since, however, chiral symmetry is an decisive feature of hadron phe- 
nomenology, it serves nevertheless to address qualitatively several questions, which 
we believe to persist to more realistic scenarios. One of these points concerns the 
question of collision and expansion time scales. We find that the expansion and 
cooling progresses relatively fast, so that the finite time, which is needed to produce 
particles, cannot be neglected. We obtain the result, that chemical and even thermal 
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equilibrium are not necessarily established during the expansion. 

Another interesting point is the production mechanism of mesons, which is stud- 
ied here for the first time. We find that mesons are created mainly within a temper- 
ature range of approximately 50 MeV below the Mott temperature. This effect can 
be traced back to a minimum of the mean free path. This feature might persist to 
a confining scenario, because it is hard to see how, even in a confining theory, pions 
could be created instantaneously. We also find that mesons are created in the bulk 
of the fireball rather than at its surface and that at late times the meson density 
follows the quark density. 

Starting from an initially spherically symmetric configuration, we observe that 
this configuration does not necessarily stay symmetric, but rather tends to form 
local density maxima at late times. All these last three statements challenge the 
applicability of hydrodynamical models. 

Using a cylindrically symmetric initial configuration, which might be more real- 
istic to encounter in an experimental situation, we find a flat rapidity distribution 
for the produced particles. The rapidity dependence of the transverse momentum is 
also flat. We observe, that in the final state mesons cover a wider range in rapidity 
than quarks, which is due to the quark mean field. 

More interesting than the SU(2) case, which is presented here, is the extension to 
577(3). There it will be possible to verify the theoretical concepts for the production 
of strangelets or strangeness destination. This will be the subject of an upcoming 
publication. 
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TABLES 



r (fm) 


y 


N u 


Nu 


N d 


N d 

a 


N n 


N n o 




K- 


3 


216 


51.5 


56.5 


52 


56 


157 


52.5 


49.5 


55 


4 


450 


106.5 


118.5 


107.5 


117.5 


422 


155 


133 


134 


5 


821 


203.5 


207 


200 


210.5 


887.5 


310 


290.5 


287 


6 


1283 


320.5 


321 


319 


322.5 


1688.5 


578 


556 


554.5 


7 


1931 


480 


485.5 


482.5 


483 


2883.5 


989 


946 


948.5 



TABLE I. Multiplicity distribution in the final state as a function of the system size. 
The columns N q and N v give the sum over all quarks and all pions, respectively. 





n q 


n n 




1 


874 


1126 


0.563 


2 


708 


1292 


0.646 


5 


538 


1462 


0.731 


10 


420 


1580 


0.790 



TABLE II. Particle multiplicity and ratio of pions to all particles for several calcula- 
tions with artificially increased hadronization cross sections. 
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FIGURES 

FIG. 1. Feynman graph for the irreducible pseudo-scalar polarization. Solid lines de- 
note constituent quarks. 

FIG. 2. Mass spectrum of the SU(2) NJL model in thermal equilibrium as a function 
of temperature. Solid line: quark mass times two, dashed line: pion mass, dot-dashed 
line: sigma mass. 

FIG. 3. Feynman graphs for the process qq — > qq. Single lines denote quarks, double 
lines mesons. The graphs for the process qq — > qq are obtained from these by a rotation 
around 90 degrees. 

FIG. 4. Feynman graphs for the process qq — > tt-k. Single lines denote quarks, double 
lines mesons. 

FIG. 5. Feynman graphs for the process tt — ► qq. Single lines denote quarks, double 
lines mesons. 

FIG. 6. Snapshot of the expansion of the fireball. Top: the system at t = 0, bottom: 
the same system at t = 25fm/c. Dark balls denote quarks, light balls denote pions. 

FIG. 7. The constituent quark mass, averaged over all solid angles, as a function of r 
at different times. 

FIG. 8. The effective temperature, averaged over all solid angles, as a function of r at 
different times. 

FIG. 9. Mean radius of the fireball as a function of time for different expansion scenar- 
ios. Full line: full interaction, dashed line: Vlasov expansion, dotted line: free expansion. 

FIG. 10. The particle density, averaged over all solid angles, as a function of r at 
different times. Solid lines: quarks, dashed lines: pions. Note the different scales on the 
vertical axes of the three rows. 
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FIG. 11. The particle multiplicity as a function of time. Solid line: quarks, dashed 
line: pions. 

FIG. 12. Time dependence of the total number of collisions per time interval. 

FIG. 13. Creation probability of pions as a function of the effective temperature. The 
Mott temperature is marked by the vertical line. 

FIG. 14. Particle density summed over quarks and mesons along the coordinate axes. 
Upper panel: t = 0, lower panel: t = 30fm/c. 

FIG. 15. Time dependence of the multiplicity for a system with cylindrically symmetric 
initial conditions. The solid line denotes quarks, the dashed line pions. 

FIG. 16. Distribution of center points in rapidity for a system with cylindrically sym- 
metric initial conditions. Upper panel: quark distribution in the initial state, lower panel: 
quark (+) and pion (x) distribution in the final state. The lines are a guide to the eye 
only. 

FIG. 17. Rapidity dependence of the mean transverse momentum for quarks in the 
initial state (+), quarks in the final state (x) and pions in the final state (*). 

FIG. 18. Energy density per particle for a fireball, in which a part of the quark mo- 
menta has been initialized according to a lower temperature. Solid line: 'hot' particles, 
dashed line: 'cold' particles. Please note the broken scale on the vertical axis. 
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